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ABSTRACT 



Joint analyses on X-ray and Sunyaev-Zel'dovich (SZ) effect of a cluster of 
galaxies can give rise to an estimate on the angular diameter distance to the 
^ ■ cluster. With the redshift information of the cluster, the Hubble constant Hq 

^ . can then be derived. Furthermore, such measurements on a sample of clusters 

with a range of redshift can potentially be used to discriminate different cosmo- 
^ . logical models. In this paper, we present statistical studies on the systematic 

^ . errors in the determination of Hq due to the triaxiality and non-isothermality of 

^ clusters of galaxies. Different from many other studies that assume artificially 

P_|. a specific distribution for the intracluster gas, we start from the triaxial model 

of dark matter halos obtained from numerical simulations. The distribution of 
the intracluster gas is then derived under the assumption of the hydrodynamic 



c§ . equilibrium. For the equation of state of the intracluster gas, both the isothermal 

and the polytropic cases are investigated. We run Monte Carlo simulations to 
^ , generate samples of clusters according to the distributions of their masses, axial 

ratios, concentration parameters, as well as line-of-sight directions. To mimic 
observations, the estimation of the Hubble constant is done by fitting X-ray and 
SZ profiles of a triaxial cluster with the isothermal and spherical /5-model. We 
find that for a sample of clusters with M = 10^^h~^ Mq and z = 0.1, the value of 
the estimated Hq is positively biased with Hq'^"'^ {estimated) ^ 1.05HQ{true) and 
Hq'"'^ [estimated) ~ 1.05Ho{true) for the isothermal case. For the polytropic case 
with 7 = 1.15, the bias is rather large with Hq'^"''^ (estimated) ~ 1.35Ho{true) 
and Hq^'^ {estimated) ~ 3HQ{true). For a mass- limited sample of clusters with 
Mii^ = lO^'-^h-^ Mq, the resuhs are similar. On the other hand, such a large 
overestimation has not been seen in real observations. It is noticed that the 
f3 value for observed clusters is within the range of 0.5 to 0.8. Considering 
only the subsample of clusters in Monte Carlo simulations with f3 in the range 
of 0.5 — 0.8, our analyses show that Hq^^ {estimated) = 1.002HQ{true) and 
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Hq'"^ (estimated) — 0.994:Ho{true) for the isothermal and polytropic cases, re- 
spectively. We further find that the value of /3 is more sensitive to the intrinsic 
asphericity of clusters of galaxies than the axial ratio of two-dimensional X-ray 
images r] is. Limiting to clusters with f3 > 0.5 essentially excludes highly aspher- 
ical clusters from the sample. From this subsample of clusters, we can get a fair 
estimate on Hq. 

Subject headings: distance scale — galaxies: cluster: general — dark matter 

1. Introduction 

The measurements of cosmological distances have played important roles in the devel- 
opment of cosmology. The discovery of the hnear distance-redshift relation for galaxies in 
the 1920s (e.g., Hubble 1929) laid the observational foundation of the big bang cosmological 
theory. The accelerating expansion of the universe revealed by the observations on Type la 
supernovae points to the existence of dark energy, which has a profound impact on cosmology 
as well as on physics (e.g., Perlmutter ct al. 1999; Riess et al. 1998). The commonly used 
distance measurements rely on various standard candles, such as Ccpheid variables. Type la 
supernovae, and the Tully-Fisher relation. The calibration of one particular standard candle 
often involves another standard candle. This is the so called distance ladder. On the other 
hand, joint analyses of X-ray and the Sunyaev-Zel'dovich effect of a cluster of galaxies can 
give rise to an estimate of the angular diameter distance to the cluster. This method is inde- 
pendent of the distance ladder, and therefore provides an important test on the consistency 
of cosmological theories. 

For clusters of galaxies, their typical mass is about 10^^ — IO^^Mq. Besides luminous 
galaxies, a large fraction of the baryonic matter in a cluster is in the form of hot intracluster 
gas. With a temperature of a few keV, the intracluster gas is fully ionized. Hot electrons emit 
strong X-ray through Bremsstrahlung processes (e.g., Rosati et al. 2002). Meanwhile, the 
electrons interact with Cosmic Microwave Background (CMB) photons and distort the CMB 
spectrum. The distortion due to the thermal motion of electrons is referred to as the thermal 
Sunyaev-Zel'dovich effect (SZ) (Sunyaev & Zel'dovich 1970, 1972). The Bremsstrahlung X- 
ray emission depends on the density of electrons through Sx cc J nlAe}i{Tc)dl, where Ue 
and Te are the number density and temperature of electrons, respectively, AgH is the X-ray 
coohng function, and the integration is along the hne of sight. For the SZ effect, we have 
ST (X J UeTedl (e.g., Carlstrom et al. 2002). From dimensional analyses, one can see that, 
up to a temperature-dependent factor, the quantity {5TY / gives rise to an estimate of the 
dimension of the cluster along the line of sight. Under the assumption of spherical symmetry. 
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this size is the same as the hnear extension of the cluster over the sky. With its measured 

angular size, one can then obtain the angular diameter distance to the cluster, and further 
the Hubble constant Hq (McHardy ct al. 1990; Birkinshaw et al. 1991; Jones et al. 1993; 
Birkinshaw & Hughes 1994; Reese et al. 2002; Reese 2004; Schmidt et al. 2004). 

Observationally, the spherical and isothermal (5 model is commonly used to describe 
the distribution of the intracluster gas. The two parameters in the model, the power index 
P and the characteristic scale Tc, are usually obtained by fitting the observed X-ray surface 
brightness with the (3 model. Thus any deviation from the sphericity and isothcrmality for 
the intracluster gas can introduce significant biases and uncertainties to the estimation of 
Hq. Studies have been performed on the effects of non- isothcrmality of the temperature 
profile, and the dumpiness and asphericity of the intracluster gas distribution (Birkinshaw 
et al. 1991; Inagaki et al. 1995; Roettiger et al. 1997; Cooray 1998; Hughes & Birkinshaw 
1998; Sulkanen 1999; Puy et al. 2000; Udomprasert et al. 2004). It is noted that the non- 
isothermality can lead to about 20 — 30% errors in the Hq determination. The dumpiness 
and asphericity may introduce about 15% errors (e.g., Inagaki et al. 1995). Most of these 
studies either assume an aspherical gas distribution of a specific form, e.g., the triaxial j3 
model (e.g., Sulkanen 1999), or concentrate on individual clusters from numerical simulations 
(e.g., Inagaki et al. 1995; Roettiger et al. 1997; Ameglio et al. 2005; Hallman et al. 2005). 
In this paper, rather than assuming an ad hoc gas distribution, we investigate the systematic 
errors starting from the triaxial mass distribution of dark matter halos. Our focus is on the 
statistical analyses for a sample of clusters generated from Monte Carlo simulations. 

From numerical simulations, Jing and Suto (2002) proposed a triaxial model with spec- 
ified statistics of the axial ratios for dark matter halos. The model has been applied to 
lensing studies (e.g., Oguri et al. 2003, 2004; Keeton et al. 2004; Tang & Fan 2005), to 
the distribution of intracluster gas (e.g., Lee & Suto 2003), and to the gas-associated X-ray 
and SZ effect of clusters of galaxies (Lee & Suto 2004; Wang & Fan 2004). In our current 
analyses on the systematic errors on Hq determination, we adopt the triaxial model for dark 
matter halos. The distribution of the intracluster gas and further the profiles of the X-ray 
surface brightness and SZ effect are derived under the assumption of hydrodynamic equi- 
librium. We consider both isothermal and polytropic equations of state for the intracluster 
gas. Regarding such distributions as 'true' distributions, the X-ray and SZ effect profiles are 
fitted with the spherical and isothermal (5 model. The two parameters (5 and Tc are obtained 
from the X-ray fitting alone as many observational analyses do. Then the estimated Hq and 
the true Hq are compared. The effects of the non-isothermality and the asphericity, as well 
as the effects of some other factors, are analyzed. To evaluate the biases and uncertainties 
statistically, we generate samples of clusters with Monte Carlo simulations. We take into 
account the distributions of the axial ratios and the concentration parameters of dark matter 
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halos. The mass function of dark matter halos is taken from Jenkins et al. (2001). Because 
of the asphericity of clusters, we also need to consider the statistics of line-of-sight directions 
in producing X-ray and SZ effect maps. 

The rest of the paper is organized as follows. In §2, we present the intracluster gas 
distribution from the triaxial model of dark matter halos. In §3, we describe the Hq de- 
termination from the joint X-ray and SZ analysis. Our results are shown in §4 and §5. §6 
contains conclusions and discussions. 



2. Density distributions and temperature profiles of the intracluster gas 

The dark matter distribution of clusters of galaxies is described by (Jing & Suto 2002) 
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where R = a{x'^ / a? + y"^ / 1)^ + / c^Y ^"^ (c < 6 < a) is the length of the major axis of a density 
contour, Rq is a characteristic scale, 5c is a dimcnsionlcss characteristic density contrast and 
Pcrit is the critical density of the universe. We take a = 1 for clusters of galaxies. 

We assume that the intracluster gas is in hydrodynamic equilibrium with the gravita- 
tional potential. Then the gas distribution is determined by 



(2) 



where Pg and Pg are the pressure and mass density of the intracluster gas, respectively, and $ 
is the gravitational potential from the total mass distribution. The specific gas distribution 
depends on its equation of state. For the isothermal gas, we have 



^ = exp[-l($-$o)], 

PgO K 



(3) 



where K — ksTg/ /imp, in which ks, Tg, /i, and nip represent the Boltzmann constant, the gas 
temperature, the mean molecular weight, and the proton mass, respectively. The subscript 
denotes the corresponding value at the central position. For the polytropic gas satisfying 
the equation of state Pg oc pj, its density and temperature distributions follow 



Pg = PgO 



Tg — TgO 
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where Kq — ksTgo/ /imp with T^o being the central temperature of the gas. 

For a cluster of galaxies, the mass fraction of the gas is about 10%. Therefore its 
gravitational potential is dominantly determined by the dark matter component. Thus 

V^* = 47rG'p . (6) 

Given the triaxial mass distribution of the dark matter in eq.(l), the gravitational potential 
can be numerically calculated. Lee and Suto (2003) derived an analytical expression of 
$ under the approximation of small asphericity for dark matter halos. The perturbative 
solution can be written in the following form 

$(«') = c[F^{u') + ^^F2{u') 

+ ^-j2 ^el{x' cos (j) — y' sin cos 9 + z' sin (p sin 6*)^ 

+el{y'sm9 + z'cos9fjF3{u')j (7) 

where u' = r'/Ro, r' = {x',y',z') with z' being the line-of-sight direction, 9 and are 
the polar coordinates of z' in the {x,y,z) coordinate system (Wang & Fan 2004), C = 
AnGScPcritRQ^, — {1 — h'^/a^y/'^ and Cc = (1 — (? ja})^!"^. Here and tc are the two 
eccentricities of the ellipsoidal dark matter halos. The three functions F\{v!\ F2{u') and 
F^{u') are given in Lee and Suto (2003). Wang and Fan (2004) further derived the X-ray 
and SZ profiles under the perturbative approximation, which are the main formulations to 
be used in our following analyses. 



3. Estimates on the Hubble constant 

In the observational analyses on X-ray and SZ effect of a cluster of galaxies, the spherical 
and isothermal (5 model is widely adopted to describe the distribution of the intracluster gas. 
With this model, the number density and the temperature of electrons are 

ne = neo[l + (r/rjr'^/', (8) 

Te = Teo, (9) 

where rieo is the central number density of electrons, Tc is the core radius, and (3 represents 
the slope of the profile. The corresponding X-ray and SZ effect profiles are 

5, = 5,o[i + ( WT^'-'^ (10) 
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and 

5T = 5T^[i + {e/e^fY'''-^^'\ (11) 

where 9 — t/Da with Da being the angular diameter distance to the cluster. In the (3 
model, the central values Sxq and 5Tq depend on ^eo, ^c, and Da- Observationally, Sxo and 
5Tq are measurable quantities, and 6c and (5 are usually obtained through the fitting to the 
X-ray profile. Thus by eliminating n^Q, the angular diameter distance Da can be estimated 
by (e.g., Reese et al. 2002) 

where G(/3) = [r(3/2/3 - l/2)/r(3/2/3)]2r(3/3)/r(3/3 - 1/2) with r(/3) being the Gamma 
function, and ^{x,Te) is the frequency-dependent factor for the SZ effect with x — hv /kBTcMB- 
The quantities /i, /c^, c, me, ctt and Tqmb are the Planck constant, the Boltzmann constant, 
the speed of light, the electron mass, the Thomson cross section, and the CMB temperature, 
respectively. The X-ray cooling function is denoted by Ae^o for = Teo, and uh = Uefie/tJ'H 
with rij = pg/fijUip for species j. From the functional form of G(/3), it is noted that eq.(12) 
can be used only with (3 > 1/3. 

For a cluster with its gas distribution being different from that of the spherical and 
isothermal f3 model, the estimated Da deviates from the true Da- Thus systematic errors 
are introduced in the Hq determination. In our analyses, we treat the gas distribution 
resulting from the triaxial model of dark matter halos as the 'true' gas distribution. For the 
isothermal case, the electron number density from the triaxial model can be approximated 
by (Wang & Fan 2004) 

Cr„.,. e^ + e? 



n. 



e 



= neoexp{--{F,{u') + ^^^F2{u') 

^— ^ ^el {x' cos (p — y' sin (p cos 9 + z' sin sin 6*) 



+ 



+el{y' sin + z' cos 0)'''^ Fs{u') } }. (13) 
The corresponding SZ effect and the X-ray profiles are 

K T 

= neQ'^{x,T,)TcMBDA{true)0o(TT ^f{0x/0o,0y/0o,Ce,eb,ec,0,(j)), (14) 



and 



" 4^(l + ^)4 -PA(^rne)go J ^nlKnd{z' / Ro) 
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= ^eo^ . , — uDA{true)9o—A^iig{e^/9o,9y/9o,Ce,eb,ec,9,(l)), (15) 
47r(l + ^j* //H 

where / = l/rieo / nf.d{z' /Rq), g — l/n^Q J nld{z'/Ro), 9q — Ro/DA{true) is the angular scale 
radius, and z is the redshift of the cluster. Then we have 

n (f ] ^ '^^^q)^ ( ^^^^ \ ^ Km^elm 1 ^(0, 0, C/K, Cb, eg, 9, 0) 

^^^'^^ 5,0 UbT,o>' 47r5?2^^^^)r^^^4(l + ^)4eom0,C^^^ ^ ^ 

Comparing DA{true) and D A{estimated) , we obtain 

DAjtrue) ^ , . 9c g{0, 0, Ce, e^, e^, ^, 0) 

D xiestimated) 9q /^(O, 0, Ce, e^, Cg, ^, 0) 

Given a cosmological model, we have Da oc therefore, 

ifo (estimated) 
ifo(true) 

^^^^o/^(O,O,Ce,e6,e„^,0)- ^ ^ 

Similarly, one can derive E2 = Ho{estimated) / Ho{true) for the polytropic case, which 
is given by 

Ho^estimated) 
HQitrue) 

6^0 /' (0, 0, Ce, 66, Cc, 9,(f))\IeO^ 

where Tew is the emission-weighted temperature, and the two functions /' and g' are defined 



as 



f 11 

f'{9J9o,9y/9o,c„eb,e,,9,<P) = j |l - — - $o]}''"c?(;27^o) (20) 
g'{9,9o,9y9o,c„eb,e,,9,<P) = J |l - - $0]} '''"'c^(^7^o) (21) 

It is seen from eq.(18) and cq.(19) that the errors on Hq estimation depend on the mass 
distribution of clusters through 9q, Cg, e^, Cc, and on the equation of state of the intracluster 
gas. Furthermore, because of the asphericity of a cluster, different observing directions also 
give rise to different results. The parameters 9c and /3 are obtained by fitting the X-ray profile 
from the triaxial model with the /3-model. Thus the values of 9c and P are dependent of the 
properties of clusters as well. Besides, because the 'true' X-ray profile for a triaxial cluster 
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does not follow the /3-profile perfectly, the fitted dc and (3 are affected by the maximum 
radius to which the fitting is applied (e.g., Navarro et al. 1995; Komatsu & Seljak 2001). It 
is noticed that Ei and E2 are very sensitive to the value of j3. Figure 1 shows the function 
of G{j3). It is seen that with (3 < 0.5, the values of Ei and E2 can be very large. For 
the polytropic case, the region used to obtain the emission weighted temperature T^w also 
matters in the estimation of Hq. 

4. Results 

In this section, we focus on the effects of the maximum fitting radius, the temperature 
gradient, and the polytropic index. We consider spherical clusters only. In the next section, 
we present our statistical analyses taking into account the asphericity of clusters. 

Theoretically, both the spherical NFW model and the triaxial model assume an infinite 
extension of the mass distribution of dark matter halos. In reality however, an astrophysical 
object always has a finite size. In our following analysis, we cut the integration range in 
eq.(14) and eq.(15) at Tint = ^r^ir. 

4.1. Effect of the maximum fitting radius 

Because the X-ray emission depends on the square of the number density of electrons, it 
decreases sharply toward the outer part of a cluster. Thus in real observations, the fitting of 
the /3-model to the X-ray profile is often performed in the inner part of a cluster where the 
signal-to-noise ratio is high enough. As we explained in the last section, different maximum 
fitting radii give rise to different values of 9c and and thus different estimates on the 
Hubble constant Hq. 

Figure 2 shows the effects of Vfu on Ei and E2 for the isothermal and polytropic cases, 
respectively. The polytropic index 7 is taken to be 1.15 in the case of E2. It is noted that 
in the isothermal case, the Hubble constant Hq is systematically overestimated with the 
maximum deviation occurring at rfu ~ O-Gr^j^- For clusters of M = 10^^h~^ Mq, the value 
of HQ{estimated) can deviate from the HQ{true) by as much as ~ 12%. The deviations are 
larger for smaller clusters. For the polytropic gas, when r/jt < 0.2r„jr, an overestimation 
is observed and the deviations arc larger for larger clusters. With rfn > 0.2rmr, the value 
of Hq is systematically underestimated with the maximum underestimation of ~ 12% at 

~ 0.6r„jj.. The smaller the clusters are, the larger the underestimations are. In the 
following analyses, we will take rfu = O.Gr^i^- For the NFW mass profile with a = 1, the 
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mass within O.Gr^jr is about 0.7 M^j^, and the average density within it is therefore about 
^Pvir ~ SOOpcrit at z = (e.g., Komatsu & Seljak 2001). Here — ^virli^'^l'^flir) ^^id 
Pcrit is the critical density of the universe. 



4.2. Effect of the temperature gradient 

When fitting the X-ray profile with the isothermal (3 model, the temperature used is the 
emission- weighted temperature T^w For the polytropic gas with a temperature distribution. 
Tew is different from the central temperature Teo, and thus the quantity {Tew/Teof^"^ in 
eq. (19) contributes an additional factor to E2. It is clearly seen that Tf.w depends on the 
temperature profile as well as on the region where the emission-weighted temperature is 
derived. 

Given the temperature profile in the polytropic gas, we can calculate the emission- 
weighted temperature through the following equation (e.g., Ascasibar et al. 2003; Rasia et 
al. 2005) 

For a spherically symmetric gas distribution, the emission-weighted temperature can be 
written as 

= j;;^ ne^T^/^r^dr ^^^^ 

Obviously, Te„, depends on the integral range r^- Figure 3 shows the influence of tt on the 
Hq estimation. In our calculations, we take 7 = 1.15, and the (3 profile fitting is done out to 
a fixed radius of rfu — O.Gr^r- when rx is varying. When rx is small, Tew/Teo is close to unity, 
and E2 > 1. With rx > 0.2rvir, E2 < 1 and reaches ~ 0.9 as rx — > r-vir. From Figure 3, it 
is seen that if we ignore the difference between and Tq, i.e., take T^w/TeQ = 1, £'2 > 1. 
Thus the underestimation appeared in the polytropic case is mainly due to Te^^/Teo < 1. In 
the analyses presented in Figure 2 as well as in the following studies, we take rx = O.Gr^, 



vir- 



4.3. The effect of the polytropic index 

Given the mass distribution of a dark matter halo, both the density and the temperature 
profiles depend on the polytropic index 7. In this subsection, we discuss how 7 affects the 
Hq estimation. 



Figure 4 presents the dependence of E2 on 7. The fitting radii for the f3 profile and the 
emission- weighted temperature are rfn = O.Gr^j^ and tt = O-Gr^ir, respectively. The mass of 
the cluster is M = 10^^h~^ M©. The result shows that E2 < 1 when 1.05 < 7 < 1.4, and 
E2 > 1 when 7 > 1.4. With 7 approaching 1, the gas becomes isothermal and E2 reaches 
~ 1.12, the result of Ei in Figure 2. At large 7 with 7 ~ 2, the value of E2 can be as high 
as about 1.6. We note that in the most theoretically probable range of 7 ~ 1.15 for the 
intracluster gas (Komatsu & Seljak 2001), the Hq is underestimated with the level of about 
10%. Our results here are consistent with those of Puy et al. (2000) and Udomprasert et al. 



Observationally, it is difficult to determine 7 precisely. The available values of 7 for 
clusters of galaxies are around 7 ~ 1.2 (e.g., Hughes et al. 1988; Finoguenov, Reiprich & 
Bohringer 2001; Pratt, Bohringer & Finoguenov 2005). Thus the polytropic index can be one 
of the factors that lead to the underestimation of Hq appeared in most of the observations. 



In this section, we discuss statistically the effects of the triaxiality of clusters on the Hq 
determination. Monte Carlo simulations are employed in the studies. 

From eq.(18) and eq.(19), it is known that Ei and E2 depend on the mass distribution 
of dark matter halos through the concentration parameter Cg, and the eccentricities e^, and 
Cc (or equivalently the axial ratios b/a and c/a), as well as on the observing direction {9, 0). 
It should be noted that the values of /5, 6c, and Tf,w/TeQ are also functions of (cg, e?,, ec, 9, cf)). 

We first analyze the effects of line-of-sight directions. Given the mass configuration of 
a cluster, the probability function of Ei can be written as 



where [^1, ^2] is the integral range of 9 within which a physically meaningful value of (j) can be 
found for a given value of Ei (Binney & de Vaucouleurs 1981). The distribution function of 
E2 has the similar form. The results are presented in Figure 5. The top and bottom panels 
are for the isothermal and polytropic cases, respectively. The position of the vertical line in 
each plot indicates the peak value of Ei[E2). Two sets of {ei),ec) are considered. They are 
(cb, Cc) = (0.2, 0.3) and (0.6, 0.8). The corresponding axial ratios are (6/a, c/a) = (0.98, 0.95) 
and (0.8, 0.6). In all the cases, the concentration parameter is taken to be the average one 
calculated from Jing and Suto (2002). The mass of the cluster is M = lO^^/i"^ M©. The 
redshift is 2; = 0.1. For the polytropic case, the value of 7 is taken to be 1.15. 



(2004) 



5. Statistical analyses on the Hq estimation in tricixial clusters 




(24) 
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With small €}, and Cc (left panels), the distribution functions of Ei and E2 are peaked 
at ~ 1.1 and ~ 0.9 with narrow extensions, respectively. For the isothermal case, the 
distribution concentrates on the side oi Ei > 1. For the polytropic case, the distribution is 
on the side of i?2 < 1- Because of the small asphericity, the results are similar to those of 
the spherical cases presented in Figure 2 with rju = O.Gr^j^. It is known that the combined 
analysis of X-ray and SZ effects of a cluster leads directly to an estimate on the line-of-sight 
dimension of the cluster. With an observed angular extension of the cluster, the angular 
diameter distance to the cluster is determined by assuming the equality of the sizes of the 
cluster parallel and perpendicular to the line of sight. Thus for a triaxial cluster oriented 
with its long axis along the line of sight, one may expect an overestimation for the angular 
diameter distance, and therefore an underestimation for Hq. The reverse trend is expected 
for a cluster with its short axis along the line of sight. Thus it seems that the distribution 
should extend to both Ei{E2) > 1 and Ei{E2) < 1, which is not seen in the left panels of 
Figure 5. The reason for this is that the above qualitative analysis considers only the effect 
of the asphericity. Taking into account the errors introduced by the differences between the 
real gas distribution and that of the isothermal (5 model, the peak values of Ei and E2 shift 
away from unity. 

With the increase of the asphericity, the distributions become broader. For the poly- 
tropic case with (cf,, Cc) = (0.6, 0.8), the values of E2 move to the side of E2 > 1. The peak 
value of E2 is around 1.8. This large overestimation is mainly due to the small (3 value 
obtained by fitting the strongly triaxial and polytropic gas distribution to the isothermal (5 
model. 

Now the distributions of Ei and E2 are investigated taking into the statistics of Cg, 
and Be- For a sample of clusters with a given mass M and at a redshift we have, for the 
isothermal case, 

(25) 

where the distributions of Cg, c/b and c/a are taken from Jing and Suto (2002) (Wang & Fan 
2004). For the polytropic case, the distribution of E2 has the similar form. 

Figure 6 shows the results for different 7. The cluster mass M = lO^^h^^ and 
z = 0.1. A clear trend is observed. As the value of 7 increases, the distribution extends 
more toward high E2. With 7 = 1.2, the tail goes well beyond £"2 = 3. The position of 
the vertical hue in each panel indicates the peak value of Ei[E2). The symbols represent 
the average values of Ei{E2). In Figure 7, we show explicitly the change of the average 
value of E2 with 7. For the isothermal case with 7 = 1, E^^^ ~ 1.06. For the case with 
7 = 1.15, £'2"^ ~ 3. It is reminded that our previous analysis on the spherical clusters with 
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7 = 1.15 gives the result of E2 ~ 0.9 (Figure 2). Thus the asphericity strongly influences 
the determination of Hq in the poly tropic case due to the small fltted value of /3. 

We present the results for clusters with different M at different z in Figure 8. The 
positions of the vertical lines represent the peak values of Ei{E2). The numbers in the 
parentheses are the average values of Ei{E2). For the isothermal case with z = 0.1, the 
average value of ^1 is 1.05, 1.06 and 1.08 for M = lO^^h-^ M©, 10^"^ /i"^ M©, and lO^^h'^ Mq, 
respectively. For M = lO^'^h^^ Mq, the values of E^^"^ at z = 0.1, 0.5 and 0.9 are respectively 
1.06, 1.08 and 1.2. In these distributions, the peak value of Ei is around 1.05. In the 
case of polytropic gas with 7 = 1.15 and z = 0.1, we have £^2^^ = 2.6,3.0 and 3.5 for 
M = 10^^h~^ Mq, 10^^/i~^ Mq, and lO^^h^^ Mq, respectively The corresponding peak 
values are 1.05, 1.35, and 1.35. With a fixed mass of M = lO^'^h-^ Mq, the peak values 
of E2 are 1.35, 1.55, and 1.55 for z — 0.1,0.5 and 0.9, respectively. The respective average 
values of E2 are 3.0,3.6 and 6.1. Thus for massive clusters at relatively high redshifts, 
the errors in Hq determination are large due to their strong asphericities. Observationally, 
however, high redshift clusters with strongly distorted X-ray/SZ images are usually not used 
in the determination of the angular diameter distance because very likely they are still in 
the merging stage. 

We further study statistically the errors for a mass-limited cluster sample with Mum = 
10^^ Mq. The probabihty function for Ei {E2) has the following form (Wang & Fan 2004) 

/ dViz) F{Ei_ 1 .1/. - ) (dn /(IM)(IM 
^^^'^ = " J{dn/dM)dMdV ' ^^^^ 

where F{Ei\M, z) is the probability density of Ei for clusters of mass M and at redshift z, 
{dn/dM)dM is the number density of clusters in the mass range of (M, M + dM), and dV 
is the volume element. The redshift range is taken to be 2; = [0, 1]. The number density of 
clusters is modeled based on the analyses of Jenkins et al. (2001). The results are shown in 
Figure 9. For the isothermal case, the distribution is peaked at £^1 ~ 1 and E'^'"'^ f« 1.06 with 
a dispersion about ±20%. For the polytropic case with 7 = 1.15, the distribution function 
has a peak at E2 ~ 1.37 and a long tail extending to high E2 values. The average value of 
E2 is £'2"^ ~ 3. Thus for a mass-hmited sample of clusters, the Hq estimation can be highly 
biased due to the combined effects of non-isothermalities and asphericities. 

There have been different studies on the systematics involved in Hq determination (e.g., 
Sulkanen 1999; Udomprasert et al. 2004). With numerical simulations, two groups recently 
analyzed possible errors in estimating Hq (Amegho et al. 2005; Hallman et al. 2005). While 
both of them noted the significant effects of the non-isothermality, they found opposite 
tendencies in the Hq determination. Ameglio et al. (2005) presented an overestimation on 
Hq, while an underestimation on Hq was reported in Hallman et al. (2005). Our analyses 
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for the spherical case show a ~ 15% underestimation in Hq with 7 in the range of 1.05 
to 1.4. The bias becomes positive for larger 7 (see Figure 4). Including the triaxiality of 
clusters, however, we find a very large bias toward overestimation in Hq with £2""^^ ~ 3 for 
7 = 1.15. The differences in the results from different groups are associated with different 
gas distributions and temperature profiles adopted in the analyses. While both Ameglio et 
al. (2005) and Hallman et al. (2005) rely on numerical simulations, we derive the properties 
of the intracluster gas from the triaxial mass distribution of dark matter halos under the 
approximation of hydrodynamic equilibrium. Even both from simulations, the obtained mass 
and temperature distributions of the intracluster gas can be quite different from one group 
to another because of the differences existed in modeling processes. It is also noted that in 
Ameglio et al. (2005), just as we have done in this paper, the parameters f3 and 6c were 
determined by fitting the f3 model to the X-ray images alone. Hallman et al. (2005), on the 
other hand, obtained the fitting parameters from X-ray and SZ images. This difference in 
the fitting technique can be an important reason to account for the different results shown by 
the two groups. The large dispersion appeared in different analyses indicates that the error 
analysis itself is model dependent. Different analyzing procedures can also lead to different 
conclusions. In our analyses, each time we change the properties of the intracluster gas, we 
obtain a new set of fitting values for /3 and 9c- Some other studies may fix the values of /3 
and 9c while changing the intrinsic gas properties (e.g., Udomprasert et al. 2004). 

Confronted with observations, the large positive bias shown in the above studies for the 
polytropic gas has not been seen in observational analyses. Most of the results on Hq derived 
from joint X-ray and SZ observations are biased low in comparison with that determined 
from other observations (e.g., Reese et al. 2002; Preedman et al. 2001), although for a couple 
of clusters the estimated values of Hq are as high as about 100 km/s/Mpc (e.g.. Mason et 
al. 2001, Jones et al. 2003). There are several issues that might be related to this apparent 
discrepancy. 

In our studies, we assume that there are no correlations between the polytropic index 
7 and the triaxiality of clusters (65, Cc). Thus, for each 7, the analyses are performed for a 
full sample of clusters with the statistics of the triaxiality given by Jing and Suto (2002). It 
is unclear, however, if 7 and {cf,, Cc) are associated with each other. In Figure 10, we plot 
the measured 7 and axial ratios 77 of the X-ray profiles for clusters of galaxies collected from 
literatures. Keeping in mind the large error bars, we see a tendency that high 7 is likely 
related with small 77, i.e., clusters with higher 7 tend to be more spherical. If the correlations 
indeed exist, the bias for polytropic and triaxial clusters can be overestimated (comparing 
Figure 7 with Figure 4). With numerical simulations, it is readily testable if there are any 
relations between 7 and the shape of the gas distribution. 
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Another issue is related to the /3 value, to which the bias is very sensitive. It is noted 
that for observed clusters, most of the f3 values are in the range of [0.5, 0.8]. Considering only 
the subsample of clusters with (3 — [0.5, 0.8], our analysis gives Ei fti 1.002 and E2 ~ 0.994. 

Observationally, clusters with strongly aspherical X-ray/SZ profiles are often avoided 
to be used in the Hq measurement. Here we study the effects of these highly elongated 
clusters on the statistical result of Hq for the sample of clusters with M = 10^^h~^ M©, 
^ = 0.1 and 7 = 1.15. In Figure 11, we show the distributions of E2 for the full sample (solid 
line), for the subsample with rj > 0.95 (dotted line) and for the subsample with (3 > 0.5 
(dash-dotted line). Here 1] is the axial ratio of the X-ray isophote contour at the virial radius 
of a cluster. It is seen that the distribution of the subsample with rj > 0.95 indeed shifts 
toward small E2 but not significantly in comparison with that of the full sample. We have 
£;f "^"(77 > 0.95) fa 1.25 and E^'"%r] > 0.95) ^ 2.5. For the full sample, the values are 
^peak ^ -j^ ^ave ^ 3 q p^j. ^^le subsamplc with P > 0.5, however, the distribution of 

E2 is narrowly around i?2 ~ 1 • We therefore immediately see that E2 depends dominantly on 
the value of f3, and the long tail in the distribution of the full sample is resulted from small 
values of f3. On the other hand, limiting rj cannot significantly affect the E2 distribution. 
The apparent reason for this is that for a given 77, a range of (3 values can exist as shown 
in Figure 12, where we apply a cut on /? > 0.33. For each rj, the maximum (3 value is 
highly correlated with the value of 77, but the scatter plot shows a large extension toward 
small p. The existence of these small P values results the long tail in the E2 distribution 
for the subsample with 77 > 0.95. We further analyze the correlations of 77 and /3 with the 
intrinsic asphericity of dark matter halos represented by and Cc. In Figure 13, we plot the 
distributions of and Cc for the three samples respectively. Clusters in the subsample with 
f3 > 0.5 are apparently biased toward spherical ones in comparison with the full sample. For 
the subsample with rj > 0.95, however, the distributions of and e^. shift only slightly to 
smaller values of and Cc- Thus the value of /? is more sensitive to e?, and than 77 is. In 
Wang and Fan (2004), we analyze the dependence of r] on the intrinsic mass distribution of 
dark matter halos. Because the distribution of the intracluster gas follows directly the profile 
of the gravitational potential, it is more spherical than the distribution of the dark matter 
is. Moreover, the shapes of X-ray/SZ images also depend on line-of-sight directions. The 
resulting 77 for the full sample of clusters discussed here is peaked around 77 ~ 0.9. In other 
words, while clusters with highly distorted X-ray/SZ images (small 7^) are strongly aspherical 
intrinsically, samples with nearly circular images (large 7/) can still contain clusters with high 
values of (e^, Cc). On the other hand, we find that the radial distribution of the intracluster 
gas is very sensitive to the intrinsic asphericity of the dark matter distribution, and thus the 
fitted P value is sensitive to (e^, Cc). Clusters with larger values of (e;,, Cc) tend to have flatter 
distributions of the intracluster gas, and therefore smaller P values. In Figure 14, we show 
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the circularized radial profiles of the X-ray surface brightness for three specific clusters all 
with T) = 0.95. Fitting the profiles with the isothermal /?-model, we get (3 — 0.4, 0.5, and 0.6, 
respectively. The respective (e?,, Cc) of underlying dark matter distributions are (0.67,0.77), 
(0.48,0.71) and (0.40,0.54). Combining the results shown in Figure 11 and Figure 13, we 
conclude that by considering intrinsically nearly spherical clusters only, we eliminate the 
large overestimates on Hq. The (3 value is a better indicator of the asphericity of underlying 
dark matter distribution than the two-dimensional axial ratio of X-ray/SZ isophote contours 
r) is. Our study presented here also opens a possibility of using the distribution of the /3 
value to probe statistically the intrinsic mass distribution of dark matter halos. 

There are other physical factors that can affect the Hq determination considerably. 
Analyses based on numerical simulations show that (e.g., Ameglio et al. 2005) excluding 
unrelaxed clusters with apparent substructures from a sample of clusters can significantly 
reduce the errors in the Hq estimation. The intrinsic anisotropics of CMB can also introduce 
a notable bias on the Hq determination (e.g., Udomprasert et al. 2004). 

6. Summctry and Discussion 

Combined analyses of X-ray and SZ effect observations provide an important means in 
determining the angular diameter distances to clusters of galaxies. In comparison with the 
method using standard candles, this determination does not depend on distance ladders, and 
thus avoids the error propagation from ladder to ladder. Furthermore, the X-ray/SZ obser- 
vations and observations on standard candles involve different systematic errors. Therefore 
important tests on the consistency of cosmological theories can be possible by comparing 
the results on the evolution of the universe from two sets of observations. Besides the de- 
termination on Hq, observations on high-redshift clusters can lead to a Hubble diagram of 
Da — z relation, which potentially can be used to constrain cosmological parameters just as 
SNe la observations do. In order to extract useful cosmological information, however, one 
needs to understand thoroughly the systematic errors involved in the method of X-ray/SZ 
observations. 

In this paper, we focus on the triaxiality and non-isothermality of clusters of galaxies. 
From the triaxial model of dark matter halos (Jing & Suto 2002), the equilibrium distribution 
of the intracluster gas is derived, which is generally aspherical and non-isothermal. The 
systematic errors in Hq estimation are evaluated with Monte Carlo simulations. 

Our main results can be summarized as follows: 

(i) Different extension adopted in fitting the /5-model to the profile of X-ray emissions 
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can introduce as much as about 12% bias in the estimated Hq value for M — 10^^h~^ M©. 
For the isothermal case, the bias is an overestimate, while it is an underestimate for the 
poly tropic gas with rfu > O-Sr^r- 

(ii) The maximum radius rr used to calculate the emission- weighted temperature in the 
polytropic case can affect Hq considerably. With rx < 0.2rvir, the value of Hq is overesti- 
mated. At larger r^, it is underestimated by about 10% with 7 = 1.15. 

(iii) In the spherical case, the effect of the non-isothermality depends on the polytropic 
index 7. With 1.05 < 7 < 1.4, Hq can be underestimated by nearly 20%. A larger 7 gives 
rise to a positive bias. As 7 is close to the adiabatic value of 1.7, the overestimate is about 
10%. 

(iv) With the increase of the asphericity of clusters, the scatters in the estimated 
Hq become large. Meanwhile, the peak value of Hq changes significantly. For the poly- 
tropic case with 7 = 1.15, the peak Hq changes from Hq^"''' {estimated) ~ 0.9HQ(true) to 
Hq'^"'^ (estimated) ^ 1.8Ho{true) as {eb,ec) increases from (0.2,0.3) to (0.6,0.8). The aver- 
age Hq^'^ (estimated) for a sample of clusters with a given mass M and at a redshift z is 
very sensitive to 7. For M = 10^^/i~^ Mq and z = 0.1, Hq^*^ (estimated) ~ 1.7Ho(true) 
at 7 = 1.05, and at 7 = 1.2, H^'"^ (estimated) ~ 3.6Ho(true). Because the asphericities of 
clusters increase with their masses and redshifts (Jing & Suto 2002), the bias on Hq increases 
with M and z as well. 

(v) For a mass-limited sample of clusters with M = 10^^h~^ M©, wc have H^"^ (estimated) 
1.06Ho(true) and the scatters are around ±20% for the isothermal case. For the polytropic 
case with 7 = 1.15, the peak value of Hq'^"''' (estimated) is about lAHo(true). There is a very 
extended tail toward high Ho(estimated) , which results an average value of Hq'"^ (estimated) ~ 
3Ho(true). 

(vi) For M — 10^^h~^ M©, z = 0.1, and 7 = 1.15, we analyze statistically the errors 
for different subsamples of clusters. For the subsample of clusters with P being in the 
range of [0.5,0.8], which is the range for most observed clusters, the average value of the 
estimated Hq is about 1.002HQ(true) and 0.994:Ho(true) for the isothermal and polytropic 
cases, respectively. We also consider a subsample of clusters with nearly circular X-ray 
images. Specifically, we select clusters with the axial ratio of two-dimensional X-ray images 
f] > 0.95. The distribution of E2 for this subsample of clusters shifts to smaller E2 side in 
comparison with the full sample. But there is still a long tail toward high values of £"2. On 
the other hand, for clusters with (3 > 0.5, the estimated E2 are around £'2 ~ 1 with ±10% 
scatters. We further show that the /3-hmited sample is clearly biased more to intrinsically 
spherical ones than the 77-limited sample is. Thus by limiting the value of /3 > 0.5, we 
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essentially select clusters that are nearly spherical intrinsically. For these clusters, we can 
obtain a fair estimate on Hq through joint X-ray/SZ analyses. 

Based on the triaxial model of dark matter halos, our analyses avoid ad hoc assumptions 
on the aspherical distribution of the intracluster gas. On the other hand, the condition of 
the hydrodynamic equilibrium used in our study is likely a simplification. The real gas 
distribution and the temperature profile of a cluster of galaxies can be more complex than 
our model can describe (e.g., Kempner et al. 2002; O'Hara et al. 2004; Reiprich et al. 2004; 
Inagaki et al. 1995; De Grandi & Molendi 2002). With dramatic advances in both X-ray 
and SZ observations, however, the knowledge on the intracluster gas will improve greatly. It 
is becoming possible observationally to describe the intracluster gas more precisely than the 
isothermal /3-model does. Therefore better determinations on Hq and even other cosmological 
parameters are highly expected from analyses on clusters of galaxies. 
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gestions. This research was supported in part by the National Science Foundation of China 
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Fig. 2. — The dependence of Ei and E2 on the maximum fitting radius Vfu. Spherical 
clusters are considered. The left panel is for the isothermal case, and the right panel is the 
polytropic case with 7 = 1.15. The pluses, asterisks and triangles are for M — 10^^h~^MQ, 
lO^^/i-^M© and 10^^ h-^M^, respectively. 



-22- 



CN 



1 .40 


_ , 1 1 1 1 , 1 1 1 , 

: 






1 


: 


1 .30 










— 


1 .20 


I \ 

= \\ 

= \\ 


— 10 










10 


\ 






1.10 


zr \\ 

= "\\ 
= \\ 


— 10 








1 .00 


z X\\ 
zr '^\\ 

I \-, \^ 










0.90 


- \" - , 

— ' ; - . _ ^ — — - 










0.80 


1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 


1 


1 





0.2 0.4 0.6 0.8 1.0 1.2 



Fig. 3. — Influence of tt on E2. Spherical clusters are considered, and the polytropic index 7 
is taken to be 1.15. The solid, dotted, and dash-dotted hues are for lO^^h'^M^, lO^^/i~^M0, 
and lO^^/i~^M0, respectively. 
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tropic cases, respectively. The vertical lines indicate the peak values of Ei{E2). 



-25- 




12 3 12 3 

E2 E2 



Fig. 6.— Probability functions F{Ei\M,z) (top left panel) and F{E2\M,z) with M = 
10^'^h~^MQ and redshift z = 0.1. In each panel, the peak value of Ei {E2) is indicated by 
the vertical line. The symbols show the average values of Ei{E2). 
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Fig. 9. — Probability functions F{Ei) and F{E2) for mass-limited cluster samples with 
-^lim — lO^^/i~^M0, The vertical lines denote the peak values of Ei and E2. The symbols 
show the average values of E-i and E2. 
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Fig. 10. — Scatter plot of 7 versus the axial ratio of the X-ray profile rj for 19 clusters. 
The triangles are from Mohr et al. (1995) (77) and Finoguenov et al. (2001) (7). The filled 
diamonds are from Mohr et al. (1995) (77) and Markevitch et al. (1999) (7). The square is 
from De Filippis et al. (2005) [f]) and Finoguenov et al. (2001) (7). The cross is from De 
Filippis et al. (2005) {r]) and Pratt et al. (2005) (7). The open diamond is from Belsole et 
al. (2005) (77 and 7). 
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Fig. 11. — The distributions of E2. The considered sample of clusters is M = 10 h Mq 
and redshift z — 0.1. The polytropic index 7 = 1.15. The solid, dotted, and dash-dotted hnes 
are for the distributions of the full sample, the subsample with rj > 0.95 and the subsample 
with P > 0.5. The positions of the vertical hnes and the symbols indicate the peak values 
and the average values of E2 of the corresponding distributions. 
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Fig. 12. — The scatter plot of f3 vs. rj for the sample of clusters with M = lO^^h 
redshift z = 0.1. The polytropic index 7 = 1.15. 
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Fig. 13. — The distributions of (left panel) and Cc (right panel) for the samples considered 
in Figure 11. 
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Fig. 14. — The circularized radial profile of X-ray emission for three specific clusters selected 
from the sample shown in Figure 12. For all the three clusters, we have r) — 0.95. The fitted 
P values are 0.4, 0.5, and 0.6, respectively. 



